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TWO ONE RANDOMIZATIONS 
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University of South Australia and Queen Mary University of London 

We investigate structure for pairs of randomizations that do not 
follow each other in a chain. These are unrandomized-inclusive, in- 
dependent, coincident or double randomizations. This involves tak- 
ing several structures that satisfy particular relations and combining 
them to form the appropriate orthogonal decomposition of the data 
space for the experiment. We show how to establish the decomposi- 
tion table giving the sources of variation, their relationships and their 
degrees of freedom, so that competing designs can be evaluated. This 
leads to recommendations for when the different types of multiple 
randomization should be used. 

1. Introduction. The purpose of this paper, and its prequel [10], is to 
establish the orthogonal decomposition of the data space for experiments 
that involve multiple randomizations [9], so that the properties of proposed 
designs can be evaluated. In [10], this was done for randomizations that 
follow each other in a chain, as in Figure 1(a). Here, analogous results to 
those in [10] are obtained for experiments in which the randomizations are 
two-to-one, as in Figure 1(b). In such randomizations, two different sets of 
objects are directly randomized to a third, as in Figures 3, 5 and 6. The 
unrandomized-inclusive, independent and coincident randomizations from 
[9] are of this type. Also covered are experiments in which the randomiza- 
tion is two-from-one in that two different sets of objects have a single set of 
objects randomized to them; that is, experiments with double randomiza- 
tions [9] [see Figure 1(c)]. 

As in [10], we always denote the set of observational units by O, so that 
the data space is the set Vh of all real vectors indexed by O. This data space 
has an orthogonal decomposition into subspaces defined by inherent factors 
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f set r 
\structure 1Z 



_T set T 
A structure Q 



(a) Two randomizations in a chain 



set fl 
structure V , 



( set r 
\structure 1Z 




I set T 
I structure Q 

(b) Two-to-one randomizations 




r set r 
V structure 1Z 



(c) Two-from-one randomizations 



Fig. 1. The three possibilities for a pair of randomizations. 



and managerial constraints. We call this decomposition the "structure" on 
fl, and identify it with the set V of mutually orthogonal idempotent matrices 
which project onto those subspaces. Thus if P G V then P is an Q x matrix, 
because its rows and columns are labelled by the elements of f2 [10]. 

In the setting of Figure 1(b), there are two other sets, T and T, which 
typically contain treatments of different types to be randomized to £1. For 
example, in Figure 3, the set of treatments (r) and the set of rootstocks (T) 
are randomized to the set of trees (O). Then Vj- is the space of all real vectors 
indexed by T, and Vr is defined similarly. Each of the sets T and T also has 
a structure defined on it, the structures being orthogonal decompositions of 
Vr and Vr, respectively. These are identified with complete sets Q and 1Z 
of mutually orthogonal idempotent matrices. 

There is an immediate technical difficulty. As first defined, a matrix Q 
in Q is not the same size as a matrix P in However, the outcome of the 
randomization of T to is a function / which allocates element f(uj) of T to 
observational unit uj. This function defines a subspace Vj- of Vq isomorphic 
to Vx . Similarly, the outcome of the randomization of T to £1 is a function g 
which allocates element g(u>) of V to observational unit uj. Thus we have a 
subspace Vf! of Vq isomorphic to Vr- From now on, we identify V-l with V?, 
and Vp with Vr- We also assume that equation (4.1) in [10] holds for both 
/ and g, so that we may regard each matrix Q in Q and each matrix R in 
1Z as an f2 x £1 matrix without losing orthogonality or idempotence. This 
condition is satisfied for all equi-replicate allocations, and for many others. 

In [10] it was seen that a standard two-tiered experiment has just two 
sets of objects, f2 and T say, typically observational units and treatments. 
To evaluate the design for such an experiment, one needs the decomposition 
of the data space Vq that takes into account both V and Q. Brien and 
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Bailey [10] introduced the notation V t> Q for the set of idempotents for 
this decomposition, and established expressions for its elements under the 
assumption that Q is structure balanced in relation to V . They exhibited 
the decomposition in decomposition tables based on sources corresponding 
to the elements of V and Q. 

For the idempotents for two sources from different tiers, such as P in 
V and Q in Q, we follow James and Wilkinson [15] in defining Q to have 
first-order balance in relation to P if there is a scalar Apq such that QPQ = 
ApqQ. If this is satisfied and Apq ^ 0, then P > Q is defined in [10] to be 
ApqPQP, which is the matrix of orthogonal projection onto ImPQ, the 
part of the source P pertaining to the source Q. The scalar Apq is called 
the efficiency factor; it lies in [0, 1] and indicates the proportion of the 
information pertaining to the source Q that is (partially) confounded with 
the source P. Furthermore, a structure Q is defined in [10] to be structure 
balanced in relation to another structure V if (i) all idempotents from Q 
have first-order balance in relation to all idempotents from V; (ii) all pairs 
of distinct elements of Q remain orthogonal when projected onto an element 
of V, that is, for all P in V and all pairs of distinct Qi and Q2 in Q, the 
product Q1PQ2 = 0. If Q is structure balanced in relation to V, and P € V , 
then the residual subspace for Q in ImP is just the orthogonal complement 
in ImP of all the spaces ImPQ: its matrix of orthogonal projection PhQ 
is given by 

(1.1) PhQ = P-^P>Q, 

where ^qgQ means summation over all Q in Q with Apq 7^ 0. 

This notation was extended in [10] to describe the decomposition for 
three-tiered experiments where the two randomizations follow each other 
in a chain, as in composed and randomized-inclusive randomizations [9] [see 
Figure 1(a)]. This involved combining the three structures V, Q and TZ 
defined on three sets of objects to yield the two equivalent decompositions 
(V \> Q) D> TZ and V > (Q t> TZ). It was seen that the idempotents of these 
decompositions could be any of the following forms: (P t> Q) [> R, P > 
(Q > R), (P > Q) h TZ, P > (Q h TZ), and P h Q, where P, Q and R are 
idempotents in V , Q, TZ, respectively. In some cases, some idempotents in 
(V D> Q) D> TZ may reduce to idempotents of the form P, P [> Q, Q, R or 
Q >R. 

In Sections 2-3 of this paper, corresponding results are obtained for 
the two-to-one randomizations: unrandomized-inclusive, independent and 
coincident randomizations. It is shown that, in addition to the decompo- 
sitions above, the following decompositions occur: V > TZ, (V t> TZ) t> Q 
and (V > Q) □ (V > TZ) , where "□" denotes "the combination of compati- 
ble decompositions" in a sense defined in Section 3. Also, the list of forms 
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of idempotents is expanded to include: P > R, P h 7Z, (P h Q) > R and 

(PhQ)h n. 

Section 4 deals with experiments having the only two-from-one random- 
ization: double randomizations. 

There are differences between different types of multiple randomization 
in the reduced forms for the above idempotents and in the efficiency factors. 
Section 5 gives recommendations for when the different types of multiple 
randomization should be used. How the results might be applied to experi- 
ments with more than three tiers is outlined in Section 6. We finish in Section 
7 with a discussion of a number of issues that arise in the decompositions 
for multitiered experiments. 

2. Unrandomized-inclusive randomizations. In an experiment with un- 
randomized-inclusive randomizations, T is randomized to £1 in an initial 
two-tiered experiment. The unrandomized-inclusive randomization involves 
a third set, T, which is randomized to 0, taking account of the result of the 
first randomization. As for randomized-inclusive randomizations, the order 
of the two randomizations is fixed. 

Two functions are required to encapsulate the results of these random- 
izations, say >T and g:f2— t-T. For oj in f2, f(oj) is the element of 
T assigned to oj by the first randomization, and g{oj) is the element of T 
assigned to oj by the second randomization. The set-up is represented dia- 
grammatically in Figure 2. 

We consider experiments in which the structure Q on T is structure bal- 
anced in relation to the structure V on Q, so that the first randomization 
gives the combined decomposition V > Q of Vq described in [10]. The sec- 
ond randomization takes account of V t> Q, both in the choice of systematic 
design and in restricting the permutations of f2 to preserve V > Q, so we 
assume that the structure 1Z on T is structure balanced in relation to V > Q- 

Put Iq = X^qgqQ' w hidi is the matrix of orthogonal projection onto 
Vy. The condition for Q to be structure balanced in relation to V can be 
written as IgPQ = ApqQ for all P in V and all Q in Q. Similarly, put 

= ^Re7e-f^i which is the matrix of orthogonal projection onto Vp- 




unrandomized in design 1 randomized in design 1 



randomized in design 2 unrandomized in design 2 unrandomized in design 2 



Fig. 2. Diagram of an experiment with two unrandomized-inclusive randomizations. 



DECOMPOSITION TABLES II. TWO-ONE RANDOMIZATIONS 



5 



Theorem 2.1. Let V , Q and 7Z be orthogonal decompositions of the 
spaces Vq, Vt and Vr, respectively, with Vy < Vq and Vr < Vq. If Q is 
structure balanced in relation to V with efficiency factors Apq, and 1Z is 
structure balanced in relation to V > Q with efficiency factors Ap^qr and 
AphQ.R, then: 

(a) 1Z is structure balanced in relation to V with efficiency matrix A-pn 
whose entries are A PR = (A P hQ,R + Z)cjeQ -Wq,r.)>' 

(b) the decomposition (V > Q) t> 7Z is 

{(P>Q)>R:PeP,QeQ,RGK,A PQ ^0, A P>Q)R + 0} 
U{(P>Q)h7e:P€P,QG Q,A PQ ^0} 
U{(PhQ)>R:PeP,RGK,A PhQiR /0} 
U{(PhQ)h7e:PG7 ? }. 

Proof, (a) Because 1Z is structure balanced in relation to the decom- 
position V \> Q, we have I^(P D> Q)R = Ap^qrR and I^(P h Q)R = 
AphQ rR, for all P in V, all Q in Q with Apq / 0, and all R in 1Z. Now, 
P = (PhQ) + E Q6Q P>Q,so 

I^PR = I W (P h Q)R + 53'l w (P > Q)R = f Api-q,r + ^'a p >q,r) R 

QeS ^ QeQ ' 

This proves that 7£ is structure balanced in relation to V with the given 
efficiency matrix. 

(b) Since 1Z is structure balanced in relation to V t> Q, we may apply the 
">" operator to elements of V > Q and 7?., to obtain 

(P > Q) > R = Apt Q iR (P > Q)R(P > Q) = A P ^ Q R (A PQ PQP)R(A PQ PQP). 

Moreover, writing ^R.e7?. ^° mean summation over RgK with Ap^qr ^ 0, 
applying equation (1.1) to P t> Q and 1Z gives 

(P > Q) h ft = P > Q - J^*(P > Q) t> R. 

Re7£ 

Similarly, 

(P h Q) > R = A p i Q R (P h Q)R(P h Q) 

and 

(P h Q) h % = P h Q - J^*(P h Q) > R. 

Thus, using Definition 4 in [10], the decomposition {V t> Q) > 1Z \s as given. 

□ 
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The expression for (V D> Q) 1> 7Z in Theorem 2.1(b) differs from that 
in equation (5.1) of [10] because (P h Q) > R is zero for composed and 
randomized-inclusive randomizations, but may not be zero for unrandomized- 
inclusive randomizations. 

For simplicity, we write the one-dimensional space for the Mean as Vo, 
with projector Po = Qo = Ro = n~ 1 3, where n = \Q\ and J is the nxn all-1 
matrix. 

As Brien and Bailey [9] show, unrandomized- inclusive randomizations are 
common in superimposed experiments. In such an experiment, it may well be 
the case that Vp n Vq~ is orthogonal to every P D> Q of the decomposition V > 
Q. In this case, the decomposition has the simpler form given by Corollary 
2.2. 

Corollary 2.2. Suppose that Q is structure balanced in relation to V 
and that 1Z is structure balanced in relation to V t> Q. If (P > Q)R = for 
all P in V, all Q in Q and all R in TZ\ {Ro}, then 

(V > Q) > K = {(P > Q) : P G V, Q G Q, A PQ + 0} 

U {(P h Q) > R : P G V, R G K, A PR ^ 0} 

U{(PhQ)hK:PGP}. 

Proof. If Api>q ; r = for all Q G Q, then Apr = Ap^R. If this is true 
for all R, then (P > Q) h ^ = P > Q. The result follows. '□ 

Lemma 2.1. Suppose that Q is structure balanced in relation to V , and 
let~P£V\ {Po}- The following conditions are equivalent. 

(i) (P > Q)Itz = for all Q in Q with A PQ / 0. 

(ii) QPR = for all Q in Q and all R in 1Z. 

(iii) I Q PI n = 0. 

If these are satsified for all P in V\ {Po} ; then Vy n Vq~ is orthogonal to 
Vr n , and all combinations of elements of T with elements of V occur 
on $1. 

Proof. If A PQ = then QP = so QPR = 0. If A PQ / then QPR = 
Apq x IgPQPI^R = Iq(P > Q)I^R. Condition (i) implies that all these 
terms are zero, which implies condition (ii). Summing QPR over all Q 
and all R gives IgPI-^,, so condition (ii) implies condition (iii). Finally, if 
Apq ^ then (P > Q)I n = A P qPQPI^ = A p 1 j PQ(IqPI7 ? .), so condition 
(iii) implies condition (i). 

Summing condition (iii) over all P in V \ {Po} gives = Ig(I-p — Po)I-r. = 
IqI-r. — IqPqI^ = (Iq — Qo)(I^ — Ro) 5 since Pq = Qo = Ro- This shows that 
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5 Viruses 



5 treatments 

\ 

10 rootstocks \ 




10 Rootstocks 



-3 Blocks 
. 10 Trees in B 



\ 

\ 



30 trees 

/ 



Fig. 3. Unrandomized-inclusive randomizations in Example 1: rootstocks are randomized 
to trees in the initial experiment; in the superimposed experiment, treatments are random- 
ized to trees taking account of the allocation of rootstocks; B denotes Blocks. 

Vy n Vq- is orthogonal to Vr H Vq . This implies that Vr D Vr = Vb, so Propo- 
sition 2 of [2] shows that the Universe is the only partition marginal to both 
T and T considered as factors on f2. Then orthogonality and Proposition 3 
of [2] show that all combinations of T and T occur on O. □ 

The conditions in Lemma 2.1 are a general form of adjusted orthogonality 
[14]. 

Example 1 (Superimposed experiment in a row-column design). The 
initial experiment in Example 10 in [9] is a randomized complete-block de- 
sign to investigate cherry rootstocks: there are three blocks of ten trees each, 
and there are ten types of rootstock. Many years later, a set of virus treat- 
ments is superimposed on this, using the extended Youden square in Table 
1. This "square" is a 3 x 10 rectangle whose rows correspond to Blocks and 
columns to Rootstocks. Each of the five treatments occurs twice in each 
Block (row), while their disposition in Rootstocks (columns) is that of a 
balanced incomplete-block design. The sets of objects for this experiment 
are trees, rootstocks and treatments. Figure 3 shows both randomizations. 

For this example, using the notation for sources in [10], but writing PMean 
as Po, the three structures are V = {Po,Pb>Pt[b]}i Q = {QchQr} and 
K = {Ro,Rv}- We have V > Q = {P > Q ,P B ,P T [B] > Qr,P T [b] ^ Q}, 
with P > Qo = Po, Pt[b] > Qr = Qr and P T[B] h Q = P T[B] - Qr. See the 
first two columns of Table 2. 

Table 1 

Extended Youden square showing the Virus Treatment for each Block- Rootstock 

combination 

Rootstocks 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Blocks I 


A 


B 


A 


C 


D 


C 


B 


E 


E 


D 


II 


D 


E 


B 


D 


E 


A 


C 


C 


A 


B 


III 


E 


A 


C 


E 


B 


D 


D 


B 


C 


A 
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The efficiency factors for the structure on treatments in relation to the 
joint decomposition of trees and rootstocks are derived from the extended 
Youden square. Viruses are orthogonal to Blocks, which means that PbRv 
0. and hence Pt[b]Rv = Rv- Viruses have first-order balance in relation 
to Rootstocks, with Ar ; v = 1/6. Hence Rv(Pt[b] t> Qr)Rv = RvQrRv = 
(l/6)Rv, and so At[b]>r,v = 1/6- Similarly, At[b]i-q,v = 5/6. Theorem 2.1 
shows that the structure on treatments is orthogonal in relation to the struc- 
ture on trees since 

^T[B],V = ^T[B]>R,V + ^T[B]hQ,V = | + t = ^ 

To obtain the full decomposition (V > Q) > 7Z, take V > Q and refine it 
by 1Z. In this experiment Vp fl Vjj 1 is not orthogonal to Vy n Vq~, because the 
Viruses source is not orthogonal to Rootstocks. This leads to nonorthogonal- 
ity between TZ and V > Q. In particular, the Viruses source is not orthogonal 
to Trees [Blocks] D> Rootstocks. Consequently, the decomposition is given by 
Theorem 2.1(b) rather than Corollary 2.2. The full decomposition of Vt rees , 
that contains six elements, one for each line in the decomposition table, is 
in Table 2: 

r (P d>Qo)>Ro,Pb, ) 

(r > q) > n = I (p T [b] > Qr) > Rv, (p T [b] > Qr) i- n, ) 

( (P T [B] I" Q) > Rv, (P T [B] I" Q) I" K ) 

with 

(P D>Q )t>Ro = Po = Qo = Ro, 
(P T[B ] > Qr) t> Rv = 6(P Tp] > Qr)R v (P T [b] > Qr) 
= 6Q r R v Qr = Qr > Rv, 
(P T [ B ] > Qr) I- K = Pt[b] > Qr - (Pt[b] > Qr) > R v 



Table 2 

Decomposition table for Example 1 



trees tier 




rootstocks tier 




treatments tier 




source 


d.f. 


source 


d.f. 


eff. 


source 


d.f. 


Mean 


1 


Mean 


1 




Mean 


1 


Blocks 


2 












Trees [Blocks] 


27 


Rootstocks 


9 


i 

6 


Viruses 
Residual 


4 

5 






Residual 


18 


5 
(i 


Viruses 
Residual 


4 
14 
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randomized unrandomized randomized 



Fig. 4. Diagram of an experiment with two independent or two coincident randomiza- 
tions. 

= Qr - Qr > Rv, 
(P t[ b] h Q) d> R V = A-f B]HS)V (P T[ B] I- Q)Rv(P T [B] I- S), 
(P T[B] HQ)hK = P T[B] HQ - (P T[B] hQ)> R V - 

As expected, this decomposition does contain a nontrivial idempotent of 
the form (P > Q) D> R. Also, unlike the chain randomizations in [10], it 
contains an idempotent of the form (P h Q) D> R. 

The efficiency factors are recorded in the decomposition in Table 2, which 
shows that the Viruses source is partly confounded with both Rootstocks and 
the part of Trees [Blocks] that is orthogonal to Rootstocks. A consequence 
of this is that four Rootstocks degrees of freedom cannot be separated from 
Virus differences. However, there are five Rootstocks degrees of freedom that 
are orthogonal to Virus differences. Further, while the Viruses source has 
first-order balance in relation to Rootstocks, the reverse is not true. 



3. Independent or coincident randomizations. For independent or coin- 
cident randomizations, two sets of objects are randomized to the third; thus 
we could have T and T randomized to Q. Two functions are needed to en- 
capsulate the results of these randomizations, say /: — > T and g: — > T. 
The set-up is represented diagrammatically in Figure 4. A particular fea- 
ture of these randomizations is that there is no intrinsic ordering of V and 
T, because neither randomization takes account of the outcome of the other. 
Associated with $7, T and T are the decompositions V, Q and 1Z. We assume 
that Q and TZ are both structure balanced in relation to V. 



r rootstocks 
t fertilizers 



r Rootstocks- 



t Fertilizers - 



C b Blocks 
r Plots in B 



t Trees in P, B 



rbt trees 



Fig. 5. Independent randomizations in Example 2: rootstocks are randomized to trees in 
such a way that all trees in each plot have a single type of rootstock; later, fertilizers are 
randomized to trees in such a way that each fertilizer is applied to one tree per plot; B, P 
denote Blocks, Plots, respectively. 
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The difference between coincident and independent randomizations is 
that, for coincident randomizations, there are sources from the two ran- 
domized tiers which are both (partly) confounded with the same source in 
the unrandomized tier. For independent randomizations this does not occur 
(apart from the Mean). 

3.1. Independent randomizations. For a pair of independent randomiza- 
tions, the two functions are randomized by two permutations chosen inde- 
pendently from the same group of permutations of f2. The precise definition 
of independence, which we were unable to give in [9], is that the conditions 
in Lemma 2.1 are satisfied, for all P in V \ {Po}, for all possible outcomes 
of the two randomizations. If ApqApr ^ then some outcomes will have 
QPR ^ 0, violating these conditions. Hence independent randomizations 
require that ApqApr = for all Q in Q and all R in 1Z unless P = Po- 
Lemma 2.1 shows that, if Q and 1Z are both structure balanced in rela- 
tion to V, then they are also structure balanced in relation to V t> 1Z and 
V D> Q, respectively, with Ap^qr = Ap^r^ = unless P = Po, Q = Qo and 
R = Ro . Therefore 

(V> Q)>K=(V>K)> Q 

= {PD>Q:Pe^,Q€Q,A PQ ^0} 
U{PhQ:PeP,PI Q ^0} 

(3.1) 

U{P>R:PGP,ReK, Apr ^ 0} 
U{PhK:PeP,PI R /0} 

u{P:PeP,PI Q = Pl7e = 0}. 

As outlined in [9], Section 8.5, wherever possible we reduce two indepen- 
dent randomizations to a single randomization. However, as noted in [9], 
Section 4.3, this is not always possible — for example, when it is not physi- 
cally possible to do them simultaneously. 

Example 2 (Superimposed experiment using split plots). Example 6 
in [9] is a superimposed experiment in which the second set of treatments 
(fertilizers) is randomized to subunits (trees) of the original experimental 
units (plots). The randomizations are independent, being carried out at 
different times and with the later one taking no account of the earlier one 
except to force fertilizers to be orthogonal to rootstocks. See Figure 5. Table 
3 shows the decomposition. 

In this example the independence of the randomizations implies that 
(Pppj > Qr) d> Rp = and so the conditions in Lemma 2.1 are satisfied. 
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Table 3 

Decomposition table for Example 2 



trees tier 


rootstocks tier 


fertilizers tier 


source 


d.f. 


source 


d.f. 


source 


d.f. 


Mean 


1 


Mean 


1 


Mean 


1 


Blocks 


6-1 










Plots [B] 


b(r - 1) 


Rootstocks 
Residual 


r - 1 
(6-l)(r-l) 






Trees [P A B] 


br(t- 1) 






Fertilizers 
Residual 


t-1 
{br-m-l) 



3.2. Coincident randomizations. For coincident randomizations, there 
are idempotents P in V \ {Po}i Q in Q \ {Qo} an d R in TZ \ {Ro} such 
that PQ and PR are both nonzero. If ImPQ and ImPR are both proper 
subspaces of ImP, then the relationship between Q and R depends on the 
choice of the two independent permutations used in randomizing T and T to 
f2; restricting one of the randomizations to preserve the relationship would 
make the multiple randomizations unrandomized inclusive rather than coin- 
cident. On the other hand, if ImPQ = ImP then ImPR is always contained 
in ImPQ. If Q is structure balanced in relation to V and ImPQ = ImP, 
then P t> Q = P and the two sources corresponding to Q and P have the 
same number of degrees of freedom. The condition for coincident random- 
izations hinted at in [9], Section 4.2, is precisely that 

, . a ^ f i n Q m 2 ano - m ^> ^ PQ an d PR are both nonzero 
then one of P \> Q and P > R is equal to P. 

A special, commonly occurring, case arises when Q and TZ can be assigned 
to the two randomized sets of objects such that the following condition is 
satisfied: 

, . for all P in V and Q in Q, if PQ and PI^ are both nonzero then 
{6 - 6) P>Q = P. 



Theorem 3.1. If Q and TZ are both structure balanced in relation to V 
and condition (3.3) is satisfied then TZ is structure balanced in relation to 
V > Q, with A p>Qi r = Apr if A PQ / and A P hQ,R = A PR i/PhQ/O. 
Moreover, the decomposition (V > Q) > TZ is 



{P > R : P e V, R G TZ, Apr + 0} 
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U{(Pr-ft):Pe7>,PI^0} 



(3.4) 



U{(P>Q):PGP,Qe Q,A P Q^O,PI^ = 0} 
U{(PhQ):PGP,PI K = 0}. 



Proof. If A PR = then PR = 0, so (P > Q)R = for all Q with A PQ + 
0, and hence (P h Q)R = 0. Suppose that PI^ / 0. Then either PI Q = 
or there is a unique Q in Q with Apq ^ 0, which satisfies P = P o Q. In 
the first case, P = Phg: therefore I^(P h Q)R = I^PR = A PR R, and 
(P h Q) > R = P > R. In the second case, I W (P > Q)R = I K PR = A PR R 
and (P > Q) t> R = P > R. □ 

Example 3 (A plant experiment). Example 5 in [9] is an experiment to 
investigate five varieties and two spray regimes. Each bench has one spray 
regime and two seedlings of each variety. See Figure 6. The sets are posi- 
tions, seedlings and regimes. The diagram includes the pseudofactor Si for 
Seedlings [Varieties], which indexes the groups of seedlings randomized to 
the different benches. Although the factor Seedlings is nested in Varieties, 
Si is not, because each of its levels is taken across all levels of Varieties. 

The Hasse diagrams displaying the structures for this experiment are in 
Figure 7. The decomposition is in Table 4, where the source Seedlings [Varieties] h 
Si is the part of Seedlings [Varieties] which is orthogonal to the source Si. 

The full decomposition of Vp OS itions in this case contains five elements and 

is 



This experiment clearly meets condition (3.3), because the only source 
for positions which is nonorthogonal to sources from both of the randomized 
tiers is the Benches source, and the five-dimensional pseudosource Si is equal 



'(P >Q )[>Ro, 
CP > Q) > K = I (Pb>Qsi)>Rr,(Pb>Qsi)I-^, 

k P P[B] > Qv, Pp[B] > Qs[V]hSi 




r 



5 Varieties 




60 positions 



60 seedlings 



10 Positions in B\ 



12 Seedlings in V 



2 regimes 



2 Regimes 



6 Benches 



Fig. 6. Coincident randomizations in Example 3: seedlings and regimes are both random- 
ized to positions; V denotes Varieties, B denotes Benches; Si and S2 are pseudofactors 
for Seedlings. 
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Universe 
1 

Benches 
6 

B A Positions 
60 



Mean 
1 
B 
5 

P[B] 
54 



Universe: Mean 
1 * 1 



Varieties: V 
5 4 




Si: Si 
6 5 



V A Seedlings: S [V] h Si 
60 50 



Universe • Mean 

Regimes • R 
2 1 



Fig. 7. Hasse diagrams for Example 3. 



to Benches. That is, Pb t> Qs x = Pb = Qs x • The other source nonorthog- 
onal to Benches is the one-dimensional source Regimes, which is a proper 
subspace of the Benches source, and so (Pb t> QsJ > Rr = Pb t> Rr = Rr- 
Consequently, the elements of the full decomposition can be written as 
follows: 



(V > Q) > Tl : 



p ,PB>RR,P B H7e, 

Pp[B] > Qv,Pp[B] > Qs[V]hSi 



On noting that P b Rr = Rr, Pp[B]Qv = Qv and P P[B ]Qs[v]hSi = Qs[V]H3n 
the decomposition further reduces to 

(V > Q) > TZ = {P , Rr, Pb - Rr, Qv, Qs[v]h Sl }■ 

Decomposition (3.4) is convenient for algorithms, because it is (V > Q) [> 
1Z, like the decompositions in Theorem 5.1(d) in [10], Theorem 2.1(b), Corol- 
lary 2.2 and equation (3.1). However, it gives the false impression that the 
decomposition of Vq must have V refined by Q, then V > Q refined by 1Z, 
suggesting that Q and 1Z have different roles. Moreover, condition (3.3) does 
not hold for all pairs of coincident randomizations. We therefore introduce 
another joint decomposition that emphasizes the symmetry between Q and 
K. 



Table 4 

Decomposition table for Example 3 



positions tier 






seedlings tier 




regimes tier 


source 


d.f. 


source 




d.f. 


source d.f. 


Mean 


1 


Mean 




1 


Mean 1 


Benches 


5 


Si 




5 


Regimes 1 



Residual 4 



Positions [Benches] 54 Varieties 4 

Seedlings [Varieties] h Si 50 



14 



C. J. BRIEN AND R. A. BAILEY 



Definition 1. Let B and C be orthogonal decompositions of the same 
space Vq- Then B is compatible with C if BC = CB for all B in B and all 
C in C. 

Lemma 2.4 in [3] shows that if B and C are compatible then the nonzero 
products BC, for B in B and C in C, give another orthogonal decomposition 
of Vh, which is a refinement of both B and C. 

Definition 2. If B and C are orthogonal decompositions of Vh which are 
compatible with each other, then the decomposition SEJC of Vh is defined 
to be 

SDC = {BC:BgS,CgC,BC^O}. 

Thus B\DC = COB. Moreover, if B and C are also both compatible with T>, 
then ,6 □ C is compatible with D, B is compatible with C □ £>, and (,6 □ C) □ D = 
BU(CUV). Hence if m are pairwise compatible then there is no 

need for parentheses in defining B1CIB2O ■ ■ ■ dB m . This decomposition 
could be referred to as u Bi combined with B2 combined with • • • combined 
with B m . v 

Lemma 3.1. //PQPRP is symmetric for all P in V , all Q in Q and 
all R in TZ, then V > Q is compatible with V > 1Z. 

Proof. If PQPRP is symmetric then PQPRP = PRPQP. Hence if 
Apq x Apr + then (P > Q)(P > R) = A p ^A pp ^PQPPRP = 

Ap^Ap R PRPPQP = (P > R)(P > Q). Thus if P > Q is defined then it 
commutes with P and with every P > R, so it commutes with P h 7Z. Sim- 
ilarly, if P > R is defined then it commutes with P h Q. Now the same 
argument shows that P h Q commutes with P h 1Z. If Pj and Pj are dif- 
ferent elements of V then Pj t> Q and Pj h Q commute with Pj > R and 
Pj h 1Z, because all products are zero. Hence V D> Q is compatible with 
V>TZ. □ 

Theorem 3.2. If the conditions in Lemma 2.1 are, or condition (3.2) 
is, satisfied, then V > Q is compatible with V [> 1Z. 

Proof. The first conditions imply that QPR = OorP = Q = R = P . 
The second implies that QPR = or PQP = A PQ P or PRP = A PR P. In 

each case, PQPRP is symmetric, so Lemma 3.1 completes the proof. □ 

Thus the decomposition {V > Q) □ {V > TZ) , which is symmetric in Q 
and TZ, can be used for coincident or independent randomizations, or for 
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unrandomized-inclusive randomizations which satisfy the conditions in Lemma 
2.1. It is the same as decomposition (3.4) for coincident randomizations when 
condition (3.3) holds, the same as the decomposition in Corollary 2.2 for 
unrandomized-inclusive randomizations when the conditions in Lemma 2.1 
hold, and the same as decomposition (3.1) for independent randomizations. 
Condition (3.2) shows that, for a pair of coincident randomizations, each 
idempotent in (V > Q) □ (V > ft) has one of the following forms: P, P t> Q, 
P > R, P h Q or P h ft. 

If a pair of coincident randomizations does not satisfy condition (3.3), then 
it may be possible to refine ft to, say, ft 2 hi such a way that ft 2 is structure 
balanced in relation to V > Q, so that the decomposition in Theorem 2.1(b) 
can be used. It is possible if R = P whenever P o R = P. 

Example 3 (Continued). As already noted, this example satisfies con- 
dition (3.3), so V \> Q is compatible with V O ft. Here 

V > Q = {P > Qo, Pb > Qsj , Pp[ B ] > Qv, Pp[ B ] > QsivjhS! } 
= {Po,Pb,Qv,Qs[v]i-Si} 

and 

V > ft = {P > R , Pb > Rr, Pb h ft, P P[B ]} 
= {P ,Rr,P b -Rr,Pp[b]}- 

Then 

(P > Q) □ CP > ft) = {Pi P b Rr, Pb(Pb - Rr), QvP P[ b], Qs[v]h Sl Pp[B]} 
= {Po,R-r,Pb - Rr, Qv, Qs[v]hSi} 
= (V > Q) > ft. 

4. Double randomizations. Double randomization is the one known type 
of two-from-one randomizations. In an experiment with double randomiza- 
tion, one set of objects is randomized to two others; thus we could have T 
randomized to T and to £1. We follow the convention that the set of obser- 
vational units is designated as ft. Two functions are needed to encapsulate 
the results of these randomizations, say /: Q — > T and g: T — > T. These two 
functions are randomized independently using two different groups of per- 
mutations. The set-up is shown in Figure 8. 

Now we obtain a subspace vJ of Vq and a subspace Vf! of Vt, both 
isomorphic to Vr- If [T[ = |T| then Vx = Vp , so we may effectively identify 
Vx, Vr and v/ ■ If |T| > |T| then we cannot identify Vy with a subspace of 
Vq without further information explicitly assigning an element of T to each 
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unrandomized randomized unrandomized 



Fig. 8. Diagram of an experiment with double randomization. 

observational unit in £1. This may not be possible (see, e.g., Figure 28 in 
[9]). Thus we shall assume that |T| = 

Associated with 0, T and T are the decompositions V, Q and 1Z. If 1Z is 
structure balanced in relation to Q and |T| = |r|, then Lemma 4.2 in [10] 
shows that Q t> 1Z = 1Z. Therefore it suffices to have 1Z structure balanced 
in relation to V . Then the overall decomposition is V \> TZ = V t> (Q > 1Z), 
which must be done from right to left. 

Example 4 (An improperly replicated rotational grazing experiment). 
Example 8 in [9] is the rotational grazing trial shown in Figure 9, with Cows 
substituted for Animals. The double randomization of Availability results in 
the assignment of Cows to Paddocks, the Cows assigned to an Availability 
forming a single herd that is used to graze all Paddocks with the same 
level of Availability. The sets of objects are observational units, paddocks 
and treatments, and the numbers of paddocks and treatments are equal, as 
required. The Hasse diagrams for treatments and observational units are like 
the middle diagram in Figure 7; that for paddocks is trivial. 

The structures on observational units, paddocks and treatments are V = 
{P ,Pc,Pr,Pc#r}, Q = {Qo,Qp} and 1Z = {R , R A , Rr, Ra#r}, respec- 
tively. This leads to the decomposition V t> (Q t> 1Z) in Table 5. It shows 
that there are no residual degrees of freedom for testing any treatment 
differences — hence the experiment being dubbed improperly replicated. 

In this case, 

Q\>1Z = K = {Q > Ro, Qp t> R A , Qp > Rr, Qp > Ra#r} 



60 observational units 



12 treatments 



f . \ 

3 Availability 

4 Rotations 

v. J 




12 paddocks 



Fig. 9. Double randomizations in Example 4'- treatments are randomized to both obser- 
vational units and paddocks. 
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Table 5 

Decomposition table for Example 4 



observational units tier 


paddocks tier 


treatments tier 




source 


d.f. 


source 


d.f. 


source 


d.f. 


Mean 


1 


Mean 


1 


Mean 


1 


Cows 


14 


Paddocks 


2 


Availability 


2 






Residual 


12 






Rotations 


3 


Paddocks 


3 


Rotations 


3 


Cows # Rotations 


42 


Paddocks 


6 


Availability # Rotations 


6 






Residual 


36 







with 

Qo > R = Ro, Q P t> R A = Ra, 

Qp > Rr = Rr, Qp > Ra#r = Ra#r- 

Also, PR is equal to either R or for all P € V and all RG?l. That is, 
1Z is orthogonal in relation to V . Therefore the complete decomposition for 
the experiment is 

V f> (O r> 7?^ = / P ° > (Q ° > Ro) ' P ° > (Qp t> Ra )' P C I" (Q » ft), 

U ' \P R C>(Qp[>RR),Pc#R>(Qpt>RA#R),Pc#R^(Q>^) 

with 

P > Qo > Ro = Po, Pc > (Qp > Ra) = Qp t> Ra = Ra, 
Pc I- (Q > n) = P c - (Qp t> R A ) = Pc - Ra, 
Pr > (Qp > Rr) = Qp > Rr = Rr, 
Pc#r > (Qp > Ra#r) = Qp > Ra#r = Ra#r, 

Pc#r \-(Q>lZ) = P c #r - (Qp > Ra#r) = Pc#r - Ra#r- 

In [9] this example was redone clS Si CcLSG of randomized-inclusive ran- 
domization, using two pseudofactors Pa and Pp for Paddocks, aliased with 
Availability and Rotations, respectively. These are required if Q itself is to 
be structure balanced in relation to V, giving a decomposition from left to 
right like the one in Section 6 in [10]. 



5. Summary. We have shown in [10] and here that, under structure bal- 
ance, the six different types of multiple randomization identified in [9] all 
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lead to orthogonal decompositions of Vq using some of the following idem- 
potents: P, P > Q, P > R, (P > Q) > R, P > (Q > R), P h Q, P h TZ, 
(P > Q) h TZ, P > (Q h 7£), (P h Q) > R and (PhQ)h 71. The differences 
between the different multiple randomizations lead to differences in the re- 
duced forms for these elements and in the efficiency factors. 

Composed randomizations. If each design is structure balanced then so 
is the composite; the decompositions V > (Q > TZ) and (V > Q) > TZ are 
equal, and so the decomposition may be done in either order; and there are 
no idempotents of the form (P h Q) [> R or (P h Q) h TZ. 

Randomized-inclusive randomizations. The structures Qi and TZi for de- 
sign 1 are refined to Q and TZ using the pseudofactors that are necessary for 
the second randomization, and then the results are the same as for composed 
randomizations. 

Unrandomized-inclusive randomizations. We must have TZ structure bal- 
anced in relation to V l> Q; use the decomposition (V > Q) E> TZ, which is 
done from left to right; if the conditions in Lemma 2.1 hold then there are 
no idempotents of the form (P > Q) l> R apart from the Mean, nor any of 
the form (P D> Q) h TZ, the decomposition V > Q is compatible with V D> TZ, 
and (V > Q) > TZ = (V > Q) □ (P > TZ). 

Independent randomizations. The conditions in Lemma 2.1 must hold; 
if both designs are structure balanced then each remains structure balanced 
after the other has been taken into account; V D> Q is compatible with V t> 
TZ; the decompositions (V > Q) > TZ, (V>TZ)> Q and (P > Q) □ (V > TZ) 
are equal; and there are no idempotents of the form (P > Q) [> R apart from 
the Mean, nor any of the form (P > Q) h TZ. 

Coincident randomizations. Condition (3.2) must hold; V 1> Q is com- 
patible with V > TZ; use the decomposition (V \> Q) \3 (V t> TZ) , whose idem- 
potents have the form P, P D> Q, P t> R, P h Q or P h TZ; if condition (3.3) 
holds, this is the same as the decomposition {V > Q) t> TZ, which is done 
from left to right; otherwise, there may be a refinement of TZ giving a left- 
to-right decomposition. 

Double randomizations. We require that |T| = [r| and that TZ be struc- 
ture balanced in relation to both Q and V , so that the decomposition is 
V >TZ = V t> (Q t> TZ), which is done from right to left. It appears that 
they can also be formulated as randomized-inclusive randomizations using 
pseudofactors to refine Q to Qi for which the left-to-right decomposition 
(V > Qi) > TZ is correct. 
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6. Structure-balanced experiments with four or more tiers. Each ex- 
periment in Sections 2-4 involves only one type of multiple randomization, 
and so involves three tiers and three structures. However, multitiered ex- 
periments are not limited to this configuration. Examples 12-14 in [9] each 
have four tiers and involve more than one type of multiple randomization. 
In general, there is the set of observational units, f2, and each randomization 
adds another set of objects with its associated tier. 

Section 7 in [10] shows how to deal with three or more randomizations 
which follow each other in a chain. Mixtures of other types of multiple 
randomization should be amenable to successive decompositions of the sort 
summarized in Section 5, so long as they are handled in the correct order. 
Thus we can use a recursive procedure in which each new structure refines 
the decomposition of Vh obtained using structures accounted for previously. 
All that is required is that each successive structure should be structure 
balanced in relation to the previous decomposition. 

One class of experiments with both two-one randomizations and chain 
randomizations consists of multiphase experiments in which different treat- 
ment factors are applied in different phases, as the following example demon- 
strates. 



Example 5 (A two-phase corn seed germination experiment). Example 
12 of [9] has the four tiers shown in Figure 10. Here we have taken the oppor- 
tunity to correct the diagram given in [9]. The 36 Lots of grain within each 
Plot should be completely randomized to Plates A Containers within each 
Interval. This will not be achieved by permuting Containers within Intervals 
and Plates within Intervals A Containers, as implied in the rightmost panel 
of Figure 10. We introduce pseudofactors Li and L2 for Lots, with nine and 
four levels, respectively, like the pseudofactors for Seedlings in Example 3. 



3 Harvesters ■ 



3 harvesters 



9 treatments 



3 Sites 
2 Blocks in S 
- 3 Plots in S, B 
36 Lots in S, B, P 

648 lots 

3 Temperatures 

3 Moistures 
v J 




18 Intervals 
4 Plates in I, C 
9 Containers in I 
648 plates 



Fig. 10. Composed and coincident randomizations in Example 5: harvesters are random- 
ized to plots; lots of grain are sampled from each plot and then randomized to plates; and 
treatments are randomized to plates; S, B, P, I, C denote Sites, Blocks, Plots, Intervals, 
Containers, respectively. 
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Table 6 

Decomposition table for Example 5 



plates tier 




lots tier 




harvesters tier 


treatments tier 


source 


d.f. 


source 


d.f. 


source 


d.f. 


source 


d.f. 


Mean 


1 


Mean 


1 


Mean 


1 


Mean 


1 


Intervals 


17 


Sites 


2 














Blocks [S] 


3 














Plots [B A S] 


12 


Harvesters 
Residual 


2 
10 






Containers [I] 


144 


Li[PAB AS] 


144 






Temperature 
Moistures 
T#M 
Residual 


2 
2 
4 
136 


Plates [C A I] 


486 


Lots[PAB AS] h 


486 











The 36 Lots must be randomly allocated to the combinations of levels of 
Li and L2, independently within each level of Sites A Blocks A Plots, so that 
neither pseudofactor corresponds to any inherent source of variation. 

At each randomization, an orthogonal design is used, so there is no diffi- 
culty in constructing the decomposition in Table 6. Here Li [P A B A S] is the 
part of the source Lots[P A B A S] which is confounded with Containers [I]. 
The source Lots[P A B A S]|- is the part of Lots[P A B A S] which is orthogonal 
to Li[P A B A S]: it is confounded with Plates [C A I]. 

Bailey [5] suggests an analysis for this example which we reproduce in the 
first three columns of Table 7(a). In this, the 3-level factors Temperature 
and Moisture have been combined into a single 9-level Treatment factor, 
the intertier interactions [9] of Sites, Harvesters and Treatments have been 
included, and the notation x is used in place of We cannot be sure, but 
it is plausible that he based this decomposition on the crossing and nesting 
relationships summarized in the formula 

(6.1) (T*H*(S/B))/Q, 

where T, H, S, B and Q represent factors for Treatments, Harvesters, Sites, 
Blocks and Plates, with 9, 3, 3, 2 and 4 levels, respectively. The sources 
derived from this are in the final column of Table 7(a), with degrees of 
freedom matching those in the preceding column. 

Revision of Table 6 along similar lines, and with pseudosources replaced 
with actual sources, yields the skeleton analysis-of- variance table in Table 
7(b). Note that, given Step 4 in Table 1 of [8], an intertier interaction will 



DECOMPOSITION TABLES II. TWO-ONE RANDOMIZATIONS 



21 



Table 7 

Skeleton analysis-of-variance tables for Example 5(a) given by Bailey [5] and (b) from 
Table 6 with intertier interactions added 



(a) 







source 


d.f. 


source from (6.1) 


Phase I: field study 




Site 


2 


S 






Experimental error (a) 


3 


B[S] 






Harvester 


2 


H 






Harvester x Site 


1 


H#S 






Experimental error (b) 


6 


H#B[S] 


Phase II: laboratory 


study 


Treatment 


8 


T 






Treatment x Site 


16 


T#S 






Experimental error (c) 


24 


T#B[S] 






Treatment x Harvester 


16 


T#H 






Treatment x Harvester x Site 


32 


T#H#S 






Experimental error (d) 


48 


T#H#B[S] 






Residual 


486 


Q[T A H A S A B] 


(b) 


plates tier 




lots tier harvesters tier 


treatments tier 


source d.f. 


source 


d.f. source 


d.f. 


source d.f. 


Mean 1 


Mean 


1 Mean 


1 


Mean 1 


Intervals 17 


Sites 


2 







Blocks [S] 3 



Plots [B AS] 12 Harvesters 2 

H#S 4 
Residual 6 



Containers[I] 144 Lots[PABAS]i 144 Treatments 8 

T#S 16 

T#H 16 

T#H#S 32 

Residual 72 



Plates[CAl] 486 Lots[PABAS] h 486 



generally occur in the right-most tier that contains a main effect in the 
interaction. Table 7(b) differs from Table 7(a) in the following ways. 

1. The rationale for the sources in Table 7(a) is unclear. We had to reverse- 
engineer it by producing formula (6.1). On the other hand, the sources 
in Table 7(b) are based on the relationships between factors within each 
tier and on the confounding between sources from different tiers. 



22 



C. J. BRIEN AND R. A. BAILEY 



2. Table 7(a) does not show, as Table 7(b) does, the successive decomposi- 
tion of the vector space indexed by the observational units. The impres- 
sion given is that there is a set of sources that arise from the field phase 
and another set that arises from the laboratory phase. 

3. Table 7(a) has four sources called "experimental error" and does not men- 
tion plates, containers, intervals, blocks, plots or lots. Hence, there is no 
indication of the sources of error variation. By contrast, each source called 
"Residual" in Table 7(b) is unambiguously identified; and the labelling 
shows that all terms are affected by variation from both phases. For ex- 
ample, the Residual for Plots [B AS], labelled Experimental error (b) in 
Table 7(a), clearly arises from variability associated with Plots within 
the Sites-Blocks combinations and variability associated with Intervals. 
Similarly, it can be seen from Table 7(b) that the Residual in Table 7(a) 
arises from variability associated with Plates and Lots. 

4. As discussed in [9], Section 7.1, the usual default is that there are no in- 
tertier interactions because such inclusions would mean that the analysis 
cannot be justified by the randomization used. It parallels the assump- 
tion of unit-treatment additivity in single-randomization experiments. 
The approach using Table 6 forces the statistician to to consult the re- 
searcher about whether intertier interactions should be included, and, if 
so, to justify them. Tables 7(a) and (b) include the intertier interactions 
of Sites, Harvesters and Treatments, which suggests that it is anticipated 
that Harvesters and Treatments will perform differently at different Sites. 

5. Even with the addition of intertier interactions, the decompositions in 
Tables 7(a) and (b) are not equivalent, and so neither are the mixed 
models underlying them. Experimental errors (c) and (d) from Table 7(a) 
are combined into the Residual with 72 degrees of freedom for Lots[P A 
B A S]i in Table 7(b). To justify an analysis based on Table 7(a), one 
would need to argue that unit-treatment interaction of Treatments with 
Blocks within Sites can be anticipated in this experiment. 

7. Discussion. 

7.1. Implications of incoherent unrandomized-inclusive randomizations. 
The phenomenon of incoherent unrandomized-inclusive randomizations is 
described in [9], Section 5.2.1. Essentially, when there has been a random- 
ization to factors that are crossed, one or more of these factors become 
nested in the second randomization. 

Consider the cherry rootstock experiment in Example 1. The trees tier 
gives an orthogonal decomposition of Vq into sources Mean, Blocks and 
Trees [Blocks] of dimensions 1, 2 and 27, respectively, in the left-hand column 
of Table 2. Similarly, the rootstocks tier decomposes Vr into sources Mean 
and Rootstocks of dimensions 1 and 9. The result of the first randomization 
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is to make the Mean sources equal and to place the Rootstocks source inside 
Trees [Blocks] , thus giving the finer decomposition of Vh shown in the middle 
column of Table 2. 

The result of the unrandomized-inclusive randomization should be to fur- 
ther decompose the decomposition resulting from the first two tiers. In the 
extended Youden square, the source Viruses is orthogonal to Blocks but par- 
tially confounded with Rootstocks, so the Viruses source defines the decom- 
position in the right-hand column in Table 2. That is, the source Viruses fur- 
ther decomposes the sources Rootstocks and the Residual for Trees [Blocks], 
as required. 

In [9] we discussed the possibility that the designer of the superimposed 
experiment ignores the inherent crossing of the factors Blocks and Root- 
stocks and randomizes Viruses to Blocks in Rootstocks in a balanced 
incomplete-block design. Then the randomizations are incoherent. The per- 
mutation group for the second randomization does not preserve the structure 
arising from the first two tiers, exhibited by the two left-most columns in Ta- 
ble 2. We can see immediately that this randomization is senseless because 
it destroys the Blocks subspace preserved by the first randomization. This 
randomization might have some appeal if no block effects had been detected 
during the 20 years of the original experiment, but then the analysis of the 
second experiment would be based on an assumed model rather than on the 
intratier structures. 

Other examples of incoherent unrandomized-inclusive randomizations are 
more complicated, and perhaps less easily detected. One is the design pro- 
posed by several authors for a split-plot experiment in which the subplot 
treatments are to be assigned using a row-column design. Example 6 illus- 
trates how consideration of the decomposition table for the proposed design 
facilitates the design process and helps the detection of incoherence. 

Example 6 (Split-plots in a row-column design). Example 11 in [9] is 
based on the design with split-plots in a row-column design given in Cochran 



4 soils 



4 Soils 



3 Benches 
-4 Plants in B 
3 Altitudes 



36 leaves i 



3 viruses 



3 Viruses 




4 Soils 

3 Benches in S 
3 Altitudes in S 



36 leaves2 



Fig. 11. Incoherent randomizations in Example 6: both soils and viruses are randomized 
to leaves, but with different structures on leaves; B denotes Benches; S denotes Soils. 
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and Cox [13], Section 7.33. Diagrams for the two randomizations are given in 
Figure 11, with leaf treatments named as viruses for clarity, soil treatments 
designated as different soils for brevity, and Altitude substituted for Layer 
so that no two factors begin with the same letter. Two diagrams are needed, 
because the assumed structure on leaves changes between the randomiza- 
tions, as shown in the two right-hand panels. At first sight, this experiment 
seems to involve unrandomized-inclusive randomizations, because soils are 
randomized to leaves in the first randomization, and then viruses are ran- 
domized to leaves, taking into account the location of the soils. However, 
the change in the assumed structure on the leaves between the two random- 
izations makes them incoherent rather than unrandomized-inclusive. 

Table 8 shows an attempt to build up a decomposition table for this 
design. The first two columns follow directly from the randomization in the 
top half of Figure 11. The third column corresponds to the leaves tier in 
the bottom half of Figure 11. When we use it to refine the decomposition 
given by the first two tiers, we find that the Soils source occurs in two tiers. 
Although this can happen in special circumstances like those in Example 
4, this is already a signal that something may be wrong. We also find that 



Table 8 

Attempted decomposition table for Example 6 



leaves i 


tier 


soils tier 


leaves2 tier* 






viruses tier 


source 


d.f. 


source d.f. 


source 


d.f. 


eff. 


source d.f. 


Mean 


1 


Mean 1 


Mean 


1 




Mean 1 


Benches 


2 




B[S] B 


2 






Plants [B] 


9 


Soils 3 


Soils 


3 










Residual 6 


B[S] h 


6 






Altitudes 


2 




A[S] A 


2 






A#B 


4 




A#B[S] a# b 


4 


? 
? 


Viruses* 2 
V#S* 2 


A#P[B] 


18 




A[S]h 


6 












A#B[S] h 


12 


? 


Viruses* 2 



? V#S* 6 
Residual 4 



'The subscripts on the sources from this tier indicate that they are the part of the source 
associated with the subscripted source in the first tier, and the subscript "h" that this is 
the part of the source orthogonal to all previous parts. 

'The partial confounding of Viruses and V#S may not have first-order balance. 
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the nesting, in this tier, of Benches within Soils gives a source B[S] with 
8 degrees of freedom. This is the sum of the previous sources Benches and 
Residual in Plants [B]; these two parts are denoted B[S]b and B[S]|- in Table 
8. Similarly, the nesting of Altitudes within Soils, in this tier, gives sources 
A[S] and A#B[S] which are each the sum of two previous sources. 

The real difficulties come when we try to incorporate the column for 
the viruses tier, because the location of the Viruses source depends on the 
outcome of the randomizations. For the outcome given in [13], Section 7.33, 
and [9], Example 11, the Viruses source does not have first-order balance in 
relation to either Altitudes # Benches or Altitudes # Plants [Benches]. The 
interaction Viruses # Soils has the same problem. 

If Altitudes # Benches is merged with Altitudes # Plants [Benches] in the 
decomposition table, then the analysis is orthogonal and is equivalent to 
that given in [13]. However, this does not allow for consistent Altitude dif- 
ferences across Plants, so it removes six spurious degrees of freedom from 
what Cochran and Cox call "Error (b)" in [13], page 310. The problem is 
that the design for the Viruses does not respect the factor relationships 
established in applying the Soils. As Yates showed in [19], if the random- 
ization respects Benches and Altitudes then a randomization-based model 
must include their interaction. 

What is needed is a design for a two-tiered experiment in which the twelve 
treatments (combinations of levels of Soils and Viruses) are randomized to 
leavesi in such a way that there is a refinement of the natural decompo- 
sition of the treatments space which is structure balanced in relation to 
Altitudes # Benches. For example, one might choose the systematic design 
in Table 9 and then randomize benches, altitudes, and plants within benches. 
In this design the twelve treatments are arranged in a (3 x 3)/4 semi-Latin 
square constructed from a pair of mutually orthogonal Latin squares of order 
3. The Viruses are arranged according to one square for soils so and s±, and 
according to the other square for S2 and S3. Theorem 5.4 in [1] shows that 
this design is the most efficient with respect to Altitudes # Benches. Let Si 



Table 9 

Proposed design for Example 6 (columns denote plants; so -S3 are different soils; 0-2 

denote viruses) 







Bench I 






Bench II 






Bench III 




Soils 


so 


81 




S3 


so 


Sl 


S2 


S3 


so 


si 


S2 


S3 


Altitude 


























Top 














1 


1 


1 


1 


2 


2 


2 


2 


Middle 


2 


2 


1 


1 








2 


2 


1 


1 








Bottom 


1 


1 


2 


2 


2 


2 














1 


1 
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be a pseudo-factor for Soils whose two levels distinguish between the first two 
and the last two levels of Soils. The design is structure balanced: Viruses 
and Viruses # Si have efficiency factor 1/2 in Altitudes # Benches, while the 
rest of the interaction Viruses $ Soils is orthogonal to Altitudes # Benches. 
It has the advantage of having 10 degrees of freedom for the Residual for 
Altitudes # Plants [Benches], two more than for the Cochran and Cox [13] 
design. 

Thus construction of the decomposition table when designing the exper- 
iment can help to detect problems with a proposed design. In this case, it 
helped to draw attention to the incoherent randomizations, to highlight the 
associated problems and to give insight into how they might be redressed. 

7.2. Other structures. All the examples in [10] and this paper are poset 
block structures, being defined by some factors and their nesting relation- 
ships, as explained in [3, 10]. More generally, a structure may be a Tjur 
structure that is defined by a family of mutually orthogonal partitions or 
generalized factors (see [17] or [2]). Again, the generalized factors are sum- 
marized in the Hasse diagram that depicts their marginality relations. There 
is one projector P for each generalized factor F, obtained from the Hasse 
diagram just as in Section 3 in [10], so that the effect of P on any vector 
is still achieved by a straightforward sequence of averaging operations and 
subtractions. It is possible for some of these projectors to be zero. Structures 
derived from tiers belong to this class. 

Another common source of structure is an association scheme [3, 7]: for 
example, the triangular scheme for all unordered pairs from a set of parental 
types, which is appropriate in a diallel experiment with no self-crosses when 
the cross is regarded as the same as the cross Then the matri- 

ces P are the minimal idempotents of the association algebra [6], and the 
corresponding subspaces are its common eigenspaces [3], Chapter 2. The 
effect of P is a linear combination of the operations of taking sums over 
associate classes. In the case of the triangular association scheme with n 
parental types, the subspaces have dimensions 1, re — 1 and n(n — 3)/2; they 
correspond to the Mean, differences between parental types and differences 
orthogonal to parental types, respectively. The decomposition 1Z% in Exam- 
ple 5 in [10] comes from an association scheme with two associate classes. 

The set of treatments in a rectangular lattice design exhibits yet another 
kind of structure [4] . Although this structure derives neither from partitions 
nor from an association scheme, the effect of each P is achieved by averaging 
and subtracting. 

The results here and in [10] apply to any structure that is an orthogonal 
decomposition of the relevant vector space, so long as each structure can 
be regarded as a decomposition of Vq. For a Tjur structure Q on a set T 
randomized to Q, condition (4.1) in [10] must hold in order for Q to be 
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regarded as an orthogonal decomposition of Vq. For structures not defined 
by partitions, it seems that we need QjX'XQj to be zero whenever ^ Qj, 
where X is the HxT design matrix. For an association scheme, this implies 
that the design must be equireplicate. The analogue of Theorem 5.1(a) in 
[10] for association schemes is given in [3], Section 7.7. 

We admit that there are relevant experimental structures, such as neigh- 
bour relations in a field or increasing quantities of dose, that are not ade- 
quately described by an orthogonal decomposition of the space. Nonetheless, 
a theory which covers designed experiments where all the structures are or- 
thogonal decompositions has wide applicability, and we limit ourselves to 
such structures here and in [10]. 

7.3. Multiphase experiments. Multiphase experiments are one of the com- 
moner types of multitiered experiment. As outlined in [9], Section 8.1, two- 
phase experiments may involve almost any of the different types of multiple 
randomizations and, as is evident from Section 5, these differ in their as- 
sumptions. 

If treatments are introduced only in the first phase, then the randomiza- 
tions form a chain, as in [10]. In [18], Wood, Williams and Speed consider a 
class of such two-phase designs for which 1Z is orthogonal in relation to the 
natural structure Q\ on the middle tier, and there is a refinement Q2 of Q\ 
such that Q2 > 1Z is structure balanced in relation to V . The results there 
are less general than ours. First, the assumptions for the second phase are in 
the nature of those for randomized-inclusive randomizations only. Second, 
the designs are restricted to those for which the design for the first phase is 
orthogonal. 

If treatments are introduced after the first phase, as in Example 5, then 
some form of two-to-one randomization is needed. Similarly, Brien and 
Demetrio [11] describe a three-phase experiment involving composed and 
coincident randomizations. 

7.4. Further work. While obtaining mixed model analyses of multitiered 
experiments has been described in [9], Section 7, and [11], it remains to 
establish their randomization analysis. The effects of intertier interactions 
on the analysis need to be investigated. We would like to establish condi- 
tions under which closed-form expressions are available for the Residual or 
Restricted Maximum Likelihood (REML) estimates of the variance compo- 
nents [16] and Estimated Generalized Least Squares (EGLS) estimates of 
the fixed effects. Also required is a derivation of the extended algorithm 
described in [12] for obtaining the ANOVA for a multitiered experiment. 

Furthermore, we have provided the basis for assessing a particular design 
for a multitiered experiment, yet general principles for designing them are 
still needed. 
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